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ABSTRACT 

We are developing the Precision Array for Probing the Epoch of Reionization (PAPER) to detect 
21cm emission from the early Universe, when the first stars and galaxies were forming. We describe the 
overall experiment strategy and architecture and summarize two PAPER deployments: a 4-antenna 
array in the low-RFI environment of Western Australia and an 8-antenna array at our prototyping site 
in Green Bank, WV. From these activities we report on system performance, including primary beam 
model verification, dependence of system gain on ambient temperature, measurements of receiver and 
overall system temperatures, and characterization of the RFI environment at each deployment site. 

We present an all-sky map synthesized between 139 MHz and 174 MHz using data from both arrays 
that reaches down to 80 mJy (4.9 K, for a beam size of 2.15e-5 steradians at 154 MHz), with a 10 
mJy (620 mK) thermal noise level that indicates what would be achievable with better foreground 
subtraction. We calculate angular power spectra {Ci) in a cold patch and determine them to be 
dominated by point sources, but with contributions from galactic synchrotron emission at lower radio 
frequencies and angular wavemodes. Although the cosmic variance of foregrounds dominates errors 
in these power spectra, we measure a thermal noise level of 310 mK at £ = 100 for a 1.46-MHz 
band centered at 164.5 MHz. This sensitivity level is approximately three orders of magnitude in 
temperature above the level of the fluctuations in 21cm emission associated with reionization. 
Subject headings: cosmology: observations, instrumentation: interferometers, radio continuum: gen- 
eral, techniques: interferometric, site testing, telescopes 



1. INTRODUCTION 

The Epoch of Reionization (EoR) marks the transi- 
tion of the primordial intergalactic medium (IGM) from 
a neutral to a highly ionized state as a result of radia- 
tion from the fir st stars and massive black holes (Loeb 
fc Barkana|2001 1. This phase transition represents a key 
benchmark in the history of cosmic structure formation 
and a major frontier of cosmic evolution yet to be ex- 
plored. Observations of Gunn-Peters on (GP) absorption 



by the IGM tow ards distant quasars ( Becker et al. |2001 
Fan et al. 12006 ) and large-scale Cosmic Microwave B, 



Jack 



ground (OlVlB) po larization from Thompson scattering 



(|Page et al.| 2007|) have constrained cosmic reionization 
redi 



between redshifts 6 < z < 14. However, deeper explo- 
ration of reionization via these probes faces fundamental 
limitations: the GP-effect saturates at low neutral frac- 
tions and CMB polarization is an integral measure of the 
Thompson optical depth to recombination. 

The most incisive probe of EoR is direct observa- 
tion of the neutra l IGM using the hydrogen 21cm line 
( [Furlanetto eral]|20d6| |Barkana fc Loeb||2005a| . The 
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rich astrophysics traced by HI and the intrinsic three- 
dimensionality o f the signal make this approach espe- 



cially appealin g (Barkana fc Loeb 2005b Loeb fc Zal- 
darriaga 2004 ) . However, the cnallenges of explor- 
mg reionization with red-shifted 21cm emission in an 
observing band below 200 MHz are daunting. EoR 
detection experiments require unprecedented levels of 
instrumental calibration and foreground characteriza- 
tion. The brightness temperatures of polarized galac- 
tic synchrotron emission, continuum point-sources, and 
galactic/extra-galactic free-free emission can exceed the 
expected ~10 mK fluctuations of the 21cm EoR signal 



by more than 5 ord ers of magnitude ( Zahn et al. 2007 
Santos et al. 20051. Wide fields-of-view (li'oVs), large 
fractional bandwidths, large numbers of antennas, sig- 
nificant RFI environments, and ionospheric variation all 
present challenges for next-generation low-frequency ar- 
rays. 

The Precision Array to Probe the Epoch of Reioniza- 
tion (PAPER) is a first-generation experiment focused on 
statistical EoR detection, but with sufficient sky cover- 
age and sensitivity to detect the very rare, largest-scale 
structures formed at the end of reionization. PAPER 
represents a focused effort to overcome the substantial 
technical challenges posed by using large, meter-wave in- 
terferometric arrays to detect the 21cm EoR signal. 

In f|2] we present the principle deployments of the PA- 
PER instrument. We then outline the architectures of 
the analog (Sj3| and digital_(Sj4l signal paths and describe 
the calibration pipeline (qSl) that is applied to the data. 
Finally, we present observational results in ^|6] that in- 
clude an all-sky map and angular power spectra mea- 
sured toward a colder patch of the synchrotron sky. 
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2. PAPER DEPLOYMENTS 

PAPER is being developed as a series of deployments 
of increasing scope to address the instrumentation, cali- 
bration, and foreground characterization challenges that 
must be surmounted in order to detect a 21cm signal 
from EoR. By characterizing and optimizing each com- 
ponent in the array with careful engineering, we hope 
to reduce the complexity of data calibration and analy- 
sis. Our staged approach allows for a systematic inves- 
tigation of observational challenges, with a capacity for 
adaptation as the characteristics of our instrument and 
of interfering foregrounds are better understood. 



TABLE 1 

PAPER Deployment Characteristics 





PGB-4 


PWA-4 


PGB-8 


PWA-64 


Deployment Date 


2004 


2007 


2006-8 


(2009) 


A^ant 




4 


4 


8 


64 


A'pol 




1 


1 


1 


4 




(str) 


.96 


.96 


0.43 


0.43 




(MHz) 


100 


150 


150 


100 


-^chan 




256 


2048 


2048 


> 2048 


Tint 




32 


7.16 


14.32 


< 7.16 


"^obs 


(days) 


1 


3 


14 


90 




(m) 


100 


150 


300 


600 



There have been three principal PAPER deployments, 
whose characteristics, along with a planned 64 antenna 
deployment, are outlined in Table [l] iVant is the num- 
ber of antennas in each deployment; TVpoi is the num- 
ber of polarization cross-multiples computed in the cor- 
relator; fJs denotes the solid angle of the primary beam 
at 150 MHz; Ai^corr is the correlated bandwidth; A^chan 
is the number of frequency channels computed over 
the correlated bandwidth; Tint is the integration time 
per visibility; Tobs is the longest continual operation; 
f^max is the approximate maximum baseline length. In 
its initial 2004 deployment at the NRAO facilities in 
Green Bank^^ the PAPER instrument consisted of four 
sleeved dipoles above planar ground-screens arranged on 
an east- west line with integrated differential amplifiers. 
This minimum-redundancy array, dubbed PGB-4, had a 
longest baseline of 100 meters and provided initial field 
experience and a means by which to evaluate our ana- 
log and digital electronics. PGB-4 established PAPER'S 
basic architecture as that of a transit array of zenith- 
pointing dipoles connected to a central correlator via 
fixed-length cables running above ground. 

The second PAPER deployment, PGB-8, began in 
2006 and underwent continual development until it was 
upgraded to PGB-16 in October 2008. PGB-8 consisted 
of eight antennas deployed on a 300m diameter circle 
(Fig. |l| at the NRAO Green Bank Galford Meadow 
site, with a small hut for rudimentary climate control 
built at the center of the array to house the receivers 
and correlator. As discussed in ^ antenna signals prop- 
agate to the hut over 75 ohm coaxial cable. The corre- 
lator itself went through several cycles of improvement 

The National Radio Astronomy Observatory (NRAO) is 
owned and operated by Associated Universities, Inc. with fund- 
ing from the National Science Foundation. 
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Fig. 1. — Illustrated above are the antenna configurations, in 
topocentric coordinates, for the PGB-8 deployment (outer, dotted) 
and the PWA-4 deployment (inner, dotted). Centered on each 
antenna are circles whose radii represent antenna elevations scaled 
by a factor of 10, with light/dark gray indicating positive/negative 
elevation with respect to a fiducial antenna. 

paralleling the development of the scalable correlator ar- 
chitecture described in S|4j In its final state, the PGB- 
8 correlator employed a packetized correlator (PaCo-8) 
that processed 2048 spectral channels across 150 MHz of 
bandwidth, computing all 4 Stokes parameters, although 
the analog system only supported one polarization. In 
early 2008, we improved the forward gains of antennas 
by adding 45° side refiectors to the ground-screens. 

Our third deployment, PWA-4, was near the Murchi- 
son Radio Observatory (MRO) site^'^ in Western Aus- 
tralia in July 2007. This path-finding array gave us first- 
hand experience with the logistics of deploying and oper- 
ating an array at a very remote site. The success of this 
deployment gave us confidence in the low RFI levels and 
data quality that can be expected from a much larger 
array at this site. PWA-4 consisted of four antennas 
with planar ground-screens (no side refiectors) arranged 
in a trapezoid pattern with a maximum baseline of 150m. 
The 4-input "Pocket Correlator" (PoCo) described in |4] 
correlated signals from these four elements. 

A fourth deployment, PGB-16, is underway. A to- 
tal of sixteen sleeved dipole antennas using ground- 
screens with side refiectors have been deployed at Gal- 
ford Meadow. Initial operations are proceeding in single- 
polarization mode using PaCo-8, but will evolve to dual- 
polarization with a PaCo-16 correlator. 

3. THE ANALOG SYSTEM 

The levels of instrumental calibration and foreground 
characterization that will be required to model and re- 
move signals interfering with an EoR detection are un- 
precedented in the 100-MHz to 200-MHz band expected 
to encompass reionization. With this in mind, we have 
taken care that each stage of our analog system (see Fig. 
[2| exhibits smooth responses as a function of frequency 

We acknowledge the Wajarri-Yamatji people of Australia as 
the Native Title Claimants of the purposed MRO lands and we 
thank them for allowing scientific activity on the site. 
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Fig. 2. — The PAPER analog signal path flows from crossed dipole elements attached to a Pseudo-Difi'erential Amplifier (PDA), through 
RG-6 coaxial cable that runs over the surface of the ground, to a Receiver Card (RC) that filters signals to a 130-MHz to 185-MHz band 
before transmitting them to the correlator. Gain elements in both the PDA and RC are sensitive to ambient temperature. An enclosure 
is under development to maintain RCs at a constant temperature, so that the only temperature-dependent gain that must be modeled is 
that of the PDA for each antenna. 

wooden framework (Fig. |3] top). This design was im- 
proved to include planar wire-mesh reflectors that attach 
to the original ground-screen (now supported by a steel 
framework) and rise outward at a 45° angle (Fig. |3] bot- 
tom) , essentially becoming a dual-polarization trough re- 
flector ( |Hall|1988| ) . This design produces a primary beam 
whose angular size more closely matches the angular size 
of colder patches in galactic synchrotron emission, ef- 
fectively reducing PAPER'S sky-noise-dominated system 
temperature. Reduction of horizon gain also mitigates 
susceptibihty to some RFI sources. We have ensured 
that these side-reflectors have not compromised the spa- 
tial and spectral smoothness of the primary beam. 

The smoothness of the primary beam as a function 
of frequency and direction enable it to be effectively 
parametrized using low-order spherical harmonics and 
frequency polynomials according to the equation: 




au{s) 



Fig. 3. — PAPER antenna elements are dual-polarization sleeved 
dipoles mounted above grounding structures. These elements have 
been designed for smooth spatial and spectral responses to facili- 
tate calibration. For the PWA-4 deployment (top) this grounding 
structure was a flat, 2mx2m ground-screen mounted on a wooden 
frame. In the later PGB-8 deployment (bottom), this structure 
was upgraded to include side reflectors that narrow the size of the 
primary beam to more closely match the size of colder patches in 
the synchrotron sky. 

and direction, thereby minimizing the number of param- 
eters needed to describe these responses and limiting the 
magnitude of errors introduced by imperfect calibration. 

3.1. Antenna Design 

PAPER antennas are designed as rugg ed dual- 

§larization versions of the sleeved dipole ( [Johnson 
33j). Crossed dipoles made from copper tubing are 
encased between two thin aluminum disks, creating a 
dual-resonance structure that broadens the antenna's fre- 
quency response. The dimensions of the tubing and disks 
have been tuned for efficient operation over a 120-MHz to 
170-MHz band. This design produces a spatially smooth 
primary beam pattern that evolves slowly with frequency. 

The antenna design includes a grounding structure 
that alleviates the gain variations that result from the 
effects of climatic conditions on the dielectric properties 
of earth ground. PGB-4 and PWA-4 deployments em- 
ployed simple wire-mesh ground planes supported by a 



^ ^ ahn{k)Yim{s) , (1) 

k=0 11=0 m=0 

where the frequency-dependence of the antenna response 
is modeled by a seventh-order polynomial and the spa- 
tial variation of each polynomial coefficient is described 
by atm coefficients of low-order spherical harmonic func- 
tions Yim{s). This flrst-order beam model smoothly in- 
terpolates to any chosen pointing and frequency, and has 
been assumed for all of the analysis presented in this pa- 
per. As shown in Figure |4] the model predicts the per- 
ceived flux densities of sources with approximately 93% 
accuracy for zenith-angles less than 45°. Errors in the 
beam model are primarily the result of the spectral re- 
sponse of the primary beam being steeper as a function 
of zenith angle than what is predicted by the model. This 
systematic effect is currently impeding the accurate mea- 
surement of spectral indices of calibrator sources. 

3.2. Pseudo-Differential Amplifier Design 

A pseudo-differential amplifier (PDA, Fig. [5]) is housed 
within the riser suspending each dipole element above its 
ground-screen. Amplifiers employ three cascaded gain 
stages where the initial stage consists of a differential 
amplifier based on NEC NE46100 bipolar junction tran- 
sistors. This circuit features an inherent direct-current 
pathway from input to ground, making the amplifier 
tolerant of electrostatic charge buil d-up that can pose 
a serious problem in dry climates (Jackson & Farrell 
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Fig. 4. — Shown in the top plot are the predicted (hncs) and 
measured (dots) responses of the PGB-8 primary beam integrated 
between 138 MHz and 174 MHz towards (top to bottom at center of 
plot) Cygnus A, Taurus A, Cassiopeia A, and Virgo A. Responses 
are relative to the primary beam's zenith response. The lower plot 
shows the residual response once a model of the primary beam has 
been subtracted. Note that the measured beam response of Taurus 
A is complicated by sidelobes of the Sun at x > 0.2. 




Fig. 5. — The dual-polarization pseudo-differential amplifier, 
pictured in its casing (upright) and with the casing removed (bot- 
tom), is housed in the riser suspending the dipole elements above 
the ground screen. The gain of this amplifier is sensitive to tem- 
perature, with a coefficient of —0.02 dB/K. 



20061. The total --32 dB gain through the PDA unit 



IS sensitive to ambient temperature with a coefRcient 
H = -0.02 dB/K over the range 5-50°C. This coeffi- 
cient enters into the voltage gain of an antenna as 



9,{t)^g.{h)H{T{t)-T{to)). 



(2) 



calibration is performed. Although not corrected for in 
the PWA-4 and PGB-8 systems, gain variation with tem- 
perature can be partially offset by recording the average 
ambient temperature of the array and applying a correc- 
tion to all recorded visibilities. We are currently devel- 
oping a system for using one antenna input as a "gain- 
o-meter" by swapping dipoles for a load attached at the 
front of the PDA unit and using the gain variation of this 
input to normalize the gains of the other antenna inputs. 
A more advanced system could also record the ambient 
temperature at each PDA to allow for per-antenna gain 
correction. 

We have estimated the receiver temperature of 
the antenna/PDA system by modeling the frequency- 
dependent auto-correlation power as a function of 
time t: 



P.{t) = \9u{t)Y 



(3) 



where gu{t) is defined in Equation pi Tsyy{t) is the tem- 
perature of a model galactic synchrotron sky weighted 
by the primary beam of an antenna; Trx is the receiver 
temperature; Ai^ is the bandwidth of a chann el, rgkv 



where T{t) is the ambient temperature as a function of 
time t and represents a fixed time at which static gain 



was modeled by scaling a 408-MHz sky map (Haslam 
|eF al. 198 2j) to 162 M Hz using a spectral index of -2.52 
( [Rogers fc Bowman 20081. As the galactic synchrotron 
sky rotates through the primary beam, the perceived Tgky 
varies, providing a modulation that enables the separa- 
tion of sky and receiver temperatures. After applying 
a first-order correction for the temperature-dependence 
of amplifier gain using measured temperatures from a 
nearby weather station, we fit an average receiver tem- 
perature of 110 K for the PGB-8 system. 

3.3. Signal Transmission 

Coaxial cables transport antenna signals from each 
PDA unit to a central processing location. Because these 
cables are not buried, they must be rugged to withstand 
harsh environmental conditions and the occasional chew 
from local fauna. We chose RG-6, 75 ohm cable with 
a polyethylene jacket for its stable propagation charac- 
teristics as a function of temperature and humidity and 
for its low cost. Signal attenuation over a 150m cable 
run is approximately 12 dB at 150 MHz, with a slope of 
-1-0.034 dB/MHz. This cable also contains a wire suitable 
for delivering DC power to PDA units. 

On the receiving end, dual-channel receiver boards con- 
sisting of amplifier stages and band-limiting filters pre- 
pare signals for digitization in the correlator. The high 
gain of these receiver cards poses a significant regenera- 
tive feedb ack concern not fou nd in heterodyne systems 
( [Slurzberg fc Oste rheld|[l96l l. This issue was mitigated 
by mounting the receivers inside a special shielded en- 
closure ( |B"radley„200"6| . The gain of the amplifiers in the 
receiver cards is also sensitive to temperature, with a co- 
efficient of —0.045 dB/K. To avoid a second temperature 
dependence in signal gain, a thermal enclosure with ther- 
moelectric heat pumps will be employed to stabilize the 
temperature of all receiver boards. The "gain-o-meter" 
approach described in §3.2| will also provide a first-order 
correction. 

4. THE DIGITAL SYSTEM 

A series of real-time digital FX correlators employ- 
ing Field-Programmable Gate Array (FPGA) processors 
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Fig. 6. — PAPER'S packetizod correlator (PaCo) follows the 
frequency— cross-multiply (FX) architecture developed by the Cen- 
ter for Astronomy Signal Processing and Electronics Research 
(CASPER). PaCo currently employs II30B boards for spectral pro- 
cessing (F engines, abo ve) and cr oss-multiplication in BEE2 boards 
(X engines, above); seelParsons et al. (2008) for details. A 10-GbE 
switch is used to route data between boards, so that data for a sub- 
set of channels from all antennas are collected at each X engine, 
where all cross-multiples are computed. 

addresses the growing digital signal processing (DSP) 
needs of visibility computation in PAPER deployments. 
These correlators are based on the architecture described 



in Parsons et al. (2008), wherein DSP engines trans- 
mit packetized data through 10-Gbit Ethernet (10-GbE) 
links to commercial switches that are responsible for 
routing data between boards. This architecture, along 
with a set of analog-to-digital converters and modular 
FPGA-based DSP hardware and a software environment 
for programming, debugging, and running them, were de- 
veloped in collaboration with the Center for Astronomy 
Signal Processing and Electronics Research (CASPER, 
Parsons et al. 2006). The flexibility of this correlator 
design shortens development time, allowing a series of 
correlators of increasing scale to be developed parallel to 
paper's incremental build-out. 

4.1. CASPER Hardware and Gateware 

The generic FX correlator architecture (Fig. |6| on 
which PAPER correlators are based consists of a set of 
modules responsible for digitizing, down-converting, and 
channelizing antenna inputs ("F" Engines), followed by 
a set of signal processing modules that cross-multiply 
all antenna and polarization samples for each frequency 
and accumulate the results ("X" Engines). The prob- 
lem of transmitting data from every F engine to every X 
engine is solved by packetizing data according to the 10- 
GbE protocol and then using a commercial switch to sort 
data streams. This approach, unique to the CASPER 
architecture for correlators of this size and bandwidth, 
avoids custom backplanes and communication protocols 
that are tailored to a single application, which are the 
traditional solution to the cross-connect problem in cor- 
relators. 

Signals from PAPER antennas are digitized by Atmel 
AT84AD001B dual 8-bit ADC chips capable of digitizing 
two streams at 1 Gsample/sec. While overrated for this 



application, this chip has been used in many CASPER 
applications and is extensively tested. The wide band- 
width of this ADC allows the entire PAPER bandwidth 
to be Nyquist-sampled as a single, real voltage and then 
digitally mixed to baseband with perfect in-phase and 
quadrature-phase components. Two ADC boards con- 
nect to each of the CASPER IBOB processing boards 
used as F engines in PAPER'S correlator. On the IBOB, 
one FPGA is responsible for the signal flow outlined in 
Figure |7]-namely, digitally down-converting four digitized 
signals, spectrally decomposing the resulting baseband 
signals using a Polyphase Filter Bank (PFB) (Crochiere 
fc Rabin^|1983 Vaidyanathan"1990p , and re-quantizmg 
spectral data to 4 bits. PAPER correlators equalize an- 
tenna spectra before re-quantization to ensure optimal 
linearity of the output power with respect to input power, 
as described in |5.1| Output data are passed briefly to 
the BEE2 boarcl (described below) where they are for- 
matted into 10-GbE packets and transmitted to a com- 
mercial switch. 

A 12-port Fujitsu CX600 10-GbE switch is employed 
for routing packets between F and X engines. Each port 
on this switch carries data for two antennas. The per- 
formance of this switch was evaluated over a period of 
16 hours, with pseudo-random data transmitted and re- 
ceived on each data link. During this time, no errors 
were detected, placing a limit on the bit-error rate of 
transmission at 2.2e-16 bits/s. Furthermore, the asyn- 
chronous correlator architecture employed by PAPER is 
tolerant of dropped packets, so errors that do occur do 
not cause system failures. 

We rely on CASPER's BEE2 boards for implementing 
X engine processing. Each of four FPGAs on a BEE2 
operates independently to unscramble packets received 
from the 10-GbE switch. Two X engine cores inside each 
BEE2 FPGA then compute cross-multiplications for all 
N antennas for 1 / iV"^ of the frequency channels and the 
results are accumulated for a selectable time interval. 
Output of accumulated data over a 100-Mbit Ethernet 
connection currently imposes a 14-second minimum in- 
tegration time. We are currently migrating this output to 
10-GbE, which will allow shorter integrations that more 
adequately resolve ionospheric fluctuations. 

4.2. PoCo and PaCo-8 

The "Pocket Correlator" (PoCo, Fig. Js] top) used 
with the PWA-4 deployment is a single IBOB system 
for correlating and accumulating four input signals (Fig. 
[?]). ADCs sample each analog input at 600 MHz and a 
75-MHz to 225-MHz band is extracted digitally. A 4- 
tap PFB decomposes this sub-band into 2048 channels, 
which are then equalized to remove per-channel scaling 
differences before being re-quantized to 4 bits. Limited 
buffer space on the IBOB permits only 1024 channels (se- 
lectable from within the 2048) to be accumulated. Ac- 
cumulated visibilities are output via serial connection to 
a host computer after being integrated for 7.16 seconds. 

The 8-antenna packetized correlator (PaCo-8, Fig. [s] 
bottom) used with the PGB-8 deployment employs four 
IBOBs and one BEE2 board that communicate through 
a Fujitsu XG700 10-GbE switch, following the architec- 
ture outlined in Figure [6] Each IBOB in this system op- 
erates identically to PoCo, but branches data from the 
equalization module to a matrix transposer for forming 
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Fig. 7. — The digital processing in an F engine of PAPER'S packetized correlator (Fig. ^ also doubles as a stand-alone "Pocket 
Correlator" (PoCo) for correlating four signal inputs. In both designs, data is re-quantized to 4 bits after spectral decomposition. To 
avoid signal distortion associated with bit truncation, signals are equalized before truncation and the inverse equalization is applied in 
post-processing. 
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Fig. 8. — For the PWA-4 deployment, a "Pocket Correlator" 
(top) was employed for correlating four single-polarization antenna 
inputs. Data were recorded in MIRIAD format on a USB disk at- 
tached to a laptop. For PGB-8 and larger deployments, a pack- 
etized correlator (bottom) connects multiple F engine processors 
through a 10-GbE switch to X engine processors and output data 
are collected by a server for recording. 

frequency-based packets. Packet data for each antenna 
are multiplexed through a point-to-point connection to 
a BEE2 FPGA and then relayed in 10-GbE format to 
the switch. A central CPU on each BEE2 board col- 
lects all visibility data after it has been accumulated for 
14.32 seconds and transmits it over 100 Mb Ethernet to 
a server where it is written to disk in MIRIAD format 



( Sault et al.||1995" l. In field deployments, correlators are 
housed inside a special shielded enclosure ( Bradley|2006 1 
to mitigate self-interference. 

5. THE CALIBRATION PIPELINE 

The greatest challenge of using low-frequency interfer- 
ometry to detect the cosmic reionization lies in character- 
izing the celestial sky and modeling the system response 
such that foregrounds to the EoR signal can be effectively 
suppressed. Various methods for estimating the effects 
of instrumental calibration on a statistical EoR detection 
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have been analyzed (Morales 2005 Morales et al. 2006 
20061, but improvmg the calibration anc 
y mstruments remains an open problem 
( [Yatawatta et al. |2008| [Bowman et al.|2007[ ). Widefields- 
ot-view, large fractional bandwidths, strong RFI envi- 
ronments, and ionospheric variation all complicate the 
calibration of next-generation low-frequency arrays such 
as PAPER, the Murchison Widefield Array (MWA), 
the Long Wavelength Array (LWA), and the LOw Fre- 
quency ARray^^ (LOFAR). Obtaining an accurate sky 
model requires accurate models of primary beams, re- 
ceiver passbands, gain variation, and array geometry. 
However, characterization of these often depends on the 
availability of an accurate sky model. The PAPER ap- 
proach to this problem has been to emphasize precision, 
with the goal of characterizing system components to 
within 1% in order to facilitate further model refinement 
via sky modeling and self-calibration. 

Although many standard tools exist for calibrating 
interferometric data, we are developing an open-source 
software project called Astronomical Interferometry in 
PYthon^'^ (AIPY) that modular izes interferometric data 
reduction to facilitate the exploration of new algorithms 
and calibration techniques. The niche targeted by this 
package are developing arrays whose evolving instrumen- 
tal characteristics require flexible, adaptable tools. AIPY 
uses the dynamically interpreted Python programming 
language to wrap together inte rfaces to packag es such 



as MIRIAD (Sault et al. 1995) and HEALPix (Gorski 



et al. 2005 1 using a common numerical array interface 



AIPY also includes object-oriented implementations of 
a wide range of algorithms use ful for radio interfer - 
ometry, includ ing W -Projection ( [Cor nwell et al."2003), 
CLEAN ( jHogbom 1974), Maximum-Entropy deconvo- 



http:/ /haystack. mit.edu/ast/arrays/mwa 
http://lwa.unm.edu 
http://www.lofar.org 
^"^ http:/ /pypi.python.org/pypi/aipy 
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Fig. 9. — Above are illustrated the fraction of usable integrations 
after fiaggin RFI for each frequency channel over 3 days of observa- 
tion with PWA-4 (light) and PGB-8 (dark). Approximately 40% 
of PGB-8 integrations are unusable owing to saturation of ADC 
inputs by strong RFI. The RFI environment at PWA is exquisite, 
with fewer channels rendered unusable by continuous narrow-band 
transmissions and almost no saturation events. 



lution (Cornwell & Evans 1985 Sault 19901, aperture 
synthesis imaging, faceted map-making, visibility sim- 
ulation, and coordinate transformations. AIPY relies 
heavily on 3rd party, public-domain modules for scien- 
tific computing, allowing its development to focus on 
interferometry-specific functionality. 

AIPY analysis centers around parametrizing interfer- 
ometric measurement equations for the instrument in 
question. However, before data from a PAPER correla- 
tor is ready to be compared to the measurement equation 
shown in Equation [4] it requires gain linearization, RFI 
excision, and crosstalk removal. 

5.1. Gain Linearization 

Mitigating the effects of data quantization in the PA- 
PER correlator requires that power levels be carefully set 
for optimal SNR and that outp ut data be converted to a 
linear power scale (see Ch. 8 of|Thompson et al.l(|2001|). 
Correction factors for gain non-linearities in digital sys- 
tems employing 1-bit and 2-bit quantization of voltage 
samples are well-known. The non-linear power response 
of the PAPER correlator, with its 4-bit correlation, is 
still signifi cant. Our approach to gain linearization is 
detailed in [Parsons et al. (|2008|. Inside the correlator, 
optimal SNK during re-quantization is ensured by apply- 
ing an equalization function that shifts data samples to 
the region of maximal linearity in the 4-bit quantization 
response curve. This equalization function is updated 
hourly to adjust for changing mean power levels as the 
galactic synchrotron sky traverses the dipole beam. The 
first steps of PAPER calibration are to deduce an average 
4-bit value from the accumulated output of the correla- 
tor, to invert the known quantiza t ion re sponse for this 



value (Fig. 13 in Parsons et al. ( |200 8[)), and then to 

sed i: 



remove the equalization function used in the correlator 

5.2. RFI Excision 

As a general solution to identifying data points that 
depart from an otherwise smooth function, we employ 
an algorithm that iteratively fits a smooth polynomial to 



visibility amplitudes, subtracts that polynomial to iden- 
tify and exclude outlying points, and repeats. This algo- 
rithm identifies outlying points by finding the mean and 
standard deviation of the logarithm of all data points and 
defining a threshold in units of the standard deviation. 
The effectiveness of this technique is improved by first 
subtracting simulated visibilities using accurate models 
of the array response and the celestial sky. 

This RFI excision algorithm is applied to spectra of 
cross-correlation data to identify narrow-band interfer- 
ence and to time-series in auto-correlation data to iden- 
tify integrations where strong RFI such as aircraft com- 
munication saturates the ADC input. Such saturation 
events result in wide-band "drop-outs" where signal is 
lost across all frequencies. In both cases, detection of an 
interference event in a given channel/integration on any 
baseline results in that data being flagged for all base- 
lines. Figure [9] illustrates the usable fraction of data av- 
eraged over a full day for the PGB-8 and PWA-4 deploy- 
ments. The RFI environment in Western Australia rep- 
resents a marked improvement over that of Green Bank, 
with fewer channels occupied by continuous RFI and al- 
most no ADC saturation events. Such events cause ap- 
proximately 40% of PGB-8 integrations to be unusable. 

5.3. Crosstalk Removal 

The inadvertent introduction of a correlation between 
two independent inputs via a cross-coupling pathway 
is unavoidable in radio interferometers. Many systems 
use phase switching to help remove crosstalk introduced 
along the signal path from an antenna to the ADC. PA- 
PER does not currently implement phase switching, al- 
though we are exploring options for doing so. However, in 
the absence of (or in addition to) phase switching, there 
are other techniques for removing crosstalk from visibil- 
ity data, provided that such crosstalk is adequately stable 
with time. PAPER'S crosstalk is modeled as an additive 
component of measured visibilities that varies slowly on 
the timescale of several hours. Because this crosstalk is 
stable on a timescale much longer than the inverse of 
the typical fringe-rate of a source over the shortest base- 
hne and because its power is concentrated at low delays, 
PAPER crosstalk is easily suppressed using fr equency- 
based time-ave rages or delay /delay-rate filters (Parsons 
fc Backer| |2009'). As with RFI, crosstalk removal is facil- 
itated by the subtraction of simulated visibilities. 

5.4. paper's Measurement Equation 

At the heart of the PAPER analysis model is a mea- 
surement equation for the single polarization visibility 
response (V,y(t)) of a baseline at channel center frequency 
V and time t: 

K(i) = G,{t) J2 ^.(5„)^.(s„)e2--(^~-"+-)+^ (4) 

n 

where s„ is a time-dependent unit vector pointing in 
the direction of a source n whose flux density is S; 
Gu{t) = gw{t)g*^{t) is the frequency-dependent baseline 
response, which is expressed in terms of the voltage gain 
of each antenna {gij,v\ Eq. ^ for constituent antennas 
indexed by i^j; A^{sn) = aii/(sn)a*jj{sn) is the direction- 
dependent response of a baseline's primary beam, ex- 
pressed in terms of the electric field response of each 
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antenna {aij,^; Eq. b represents the baseline vector 
separating antenna positions in units of light propagation 
time; t is the relative electronic delay of antenna signals 
into the correlator; and <f)\s a, relative instrumental phase 
between antenna signals. 

The design of our antenna and analog electronics leads 
to antenna-based gains {aij,gij) that can be modeled by 
smooth functions in time and frequency with a modest 
number of parameters. The degrees of freedom in these 
parameters are over-constrained by wide-band visibility 
data, even when using a sky model that includes only a 
small number of point sources. However, the larger FoV 
of PAPER elements has complicated early calibration by 
decreasing the extent to which a single source dominates 
the correlated signal between antennas. Without isola- 
tion of sources, self-calibration cannot proceed as a direct 
computation using raw data, but rather must proceed by 
fitting models of the array and the observed sky to re- 
move baseline-dep endent interference patterns ( Cornwell 
& Fomalont 1989 1. To address this problem, we employ a 
two-tiered calibration process wherein initial coarse cal- 
ibration that employs source-isolation techniques is fol- 
lowed by multi-source, least-squares fitting to visibilities 
for the parametrization described in Q. 

5.5. First-Order, Time-Independent Self- Calibration 

Phase calibration involves solving for the position and 
electronic signal delay associated with each antenna. The 
parameter space associated with phase calibration con- 
tains many local minima separated by wavelength incre- 
ments projected toward calibrator sources. To speed the 
phase calibration process, we rely on theodolite-surveyed 
antenna positions. When initial antenna positions in- 
clude errors greater than a wavelength, we have devel- 
oped a procedure for fitting positions that is reason- 
ably robust against converging to non-global minima. 
First, a single source is isolated in data either by se- 
lecting a time when one source is d ominant, or by us- 
ing c oarse delay /delay-rate filtering (Parsons & Backer 



2009|). Baseline components bx,by and a throw-away 



phase term are then fit for a pair of nearby antennas 
using a narrow set of frequency channels and a range 
of time over which fringes wrap several times. Next, 
data including a few other sources are used to separate 
b^ out of the remaining phase term. As a final step in 
single-baseline phase calibration, wide-bandwidth data 
are used to account for the final phase as a combina- 
tion of electronic signal delay (t) and a small, frequency- 
independent phase term (</)). 

Following accurate phase calibration for a subset of 
antennas, another antenna is added and the procedure 
outlined above is repeated using data for all baselines 
that connect this antenna to the calibrated subset. With 
phase calibration for four antennas, preliminary antenna- 
based gain calibration is possible. PAPER visibilities 
have been flux-calibrated using Cygnus A with a value 
of 1.09e4 Jy at 150 MHz and a spe ctral index of —0.69, 



which are derived from the data in Baars et al. (19771 



The beam model described in p.l| provided a tirst-order 



correction for the effects of the primary beam on the 
Cygnus A spectrum. In addition to an initial gain model 
that incorporated the measured responses of the analog 
and digital filters in the PAPER pipeline, fourth-order 
polynomials in frequency were fit for the bandpass func- 



tions (gii/) of antennas in the array. 

At each stage, the residuals remaining after flux- 
calibration can help to highlight baselines with poorer flt 
scores-a hallmark of convergence to a non-global mini- 
mum. Once all problematic baselines have been elimi- 
nated, simultaneous fitting of phase and gain parame- 
ters can proceed until time-variable effects such as iono- 
spheric distortion and changing primary beam response 
towards strong sources dominate residual visibilities. 

5.6. Excision of Strong Point-Sources 

Strong point-sources, while useful for initial calibra- 
tion, complicate further calibration by degrading sensi- 
tivity toward other positions on the sky. This is espe- 
cially true for small arrays, since the sidelobe response 
toward a source scales approximately as ^/Nvis, where 
Nyig is the number of independent locations in the uv- 
plane that are being phased and summed. For PWA-4 
and PGB-8, sources stronger than a few hundred Jy can 
dominate the response of a synthesized beam, even when 
using earth rotation synthesis. 

Deep removal of the strongest sources is essential for 
improving primary beam models using sources at a va- 
riety of declinations. A time-invariant model of source 
position, flux density, and spectral index is moderately 
effective at removing these sources, depending on the ac- 
curacy of current beam models and the activity of the 
ionosphere. For PGB-8, we were able to achieve ~95% 
suppression of sources using static modeling. However, 
the residual flux densities of the strongest flve sources 
(Sun, Cygnus A, Cassiopeia A, Virgo A, Taurus A) were 
still strong enough to introduce noticeable imaging side- 
lobes. Improving source removal required fitting for off- 
sets in source position and flux spectrum as a function of 
time. We attribute these variations to ionospheric refrac- 
tion (Fig. 10 1 and defects in the primary beam model, 
respectively. All of these components were necessary to 
achieve source suppression at a level of le-3 to 2e-4, de- 
pending on a source's position within the primary beam. 

6. OBSERVATIONAL RESULTS 

Results were obtained using 3-day observations with 
a single east-west linear polarization for both PWA-4 
and PGB-8 deployments. For both sets of data, the five 
strongest sources visible from the northern hemisphere 
(Sun, Cygnus A, Cassiopeia A, Virgo A, and Taurus 
A) w ere removed following the procedure described in 
j |5.6[ Residuals associated with these sources refiect de- 
viations of their spectra (as viewed through an imper- 
fectly calibrated primary beam) from a strict power-law. 
The Sun was further suppressed via the application of 
a delay/fringe-rate filter (j |5.5| . Although these filters 
potentially corrupt wider areas of the map, the strength 
and variability of the Sun is such that their application 
greatly improves imaging of the daytime sky. 

6.1. PWA-4/PGB-8 All-Sky Map 

The all-sky map in Figure [13] illustrates average fiux 
density per beam area over a band between 138.8 MHz 
and 174.0 MHz. The northern hemisphere was imaged 
using PGB-8 data with peak beam response at declina- 
tion +38.5°. The southern hemisphere was imaged with 
data from PWA-4, which employed both a more compact 
array configuration and a broader primary beam. Peak 
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Fig. 10. — Illustrated above are the magnitudes of the angular offsets, primarily in the zenith direction, of Cygnus A (gray dots) and 
Cassiopeia A (black dots) fit as a function of time relative to their cataloged positions for 3 days of observations (top to bottom) with 
PGB-8. Both sources exhibit greater variability in position v^hen the Sun wfas below the horizon (gray shaded area), presumably as a 
result of greater ionospheric variability at those times. The offsets of both sources also exhibit a smoothly v arying component that repeats 
day-to-day, which may be attributed to refraction from spherical and wedge components of the ionosphere | |Komesaroff|1960| l. 



beam response in the southern hemisphere hes at dec- 
hnation —26.7°. Imaging over the full 35.2-MHz band- 
width was performed by generating full maps in 1.46- 
MHz intervals and then summing maps. Each map was 
generated using 200 facets, with phase centers equally 
spaced around the sp here. While each ima ged facet em- 
ployed W-projection ( Cornwell et al.[20"03 l to correct for 
the curvature of the sky, facet imaging was nonetheless 
required to incorporate differences in data weighting be- 
tween phase centers. This weighting optimizes SNR by 
accounting for the array's direction-dependent gain to- 
ward that location as it drifts through the primary beam. 

Thermal noise and point-source sidelobes both con- 
tribute to the noise level seen in this all-sky map. As 
mentioned in S|6] only the strongest five sources have 
been removed from this map prior to imaging. Image- 
domain CLEAN deconvolution was used to achieve mod- 
est suppression of sidelobes from other sources. The 
complexities of wide-field imaging with drift-scan data 
make further post-imaging deconvolution impractical, 
since CLEAN components of gridded visibilities may not 
be directly applied to ungridded visibility data. Our plan 
is to use imaging as a pathway for identifying sources that 
may then be fit and removed from raw visibilities. With 
this in mind, we have devoted little additional effort to- 
ward suppressing sidelobes of imaged sources. From the 
map in Figure 13 we measure an RMS flux density of 
noise between point sources of 80 mJy (4.9 K using a syn- 
thesized beam are a of 2 .15e-5 steradians at 156.4 MHz). 
As we discuss in 96.21 and in 36.31 this noise level ex- 



ceeds thermal noise by a factor of eight, and is dominated 
by point-source sidelobes, with only a small contribution 
from galactic synchrotron emission at lower frequencies. 

The position-dependent response of primary beams in 
each array is evidenced by changing noise levels in the 
map as a function of declination. Noise levels in the map 
also change with right-ascension because system temper- 
ature depends on the average temperature of galactic 
synchrotron emission across the primary beam. To il- 
lustrate the levels of thermal noise in the all- sky map, 
we repeated the map-making process detailed above, but 
with alternating integrations added with opposite signs, 
so that celestial sources with slow time variation were 



heavily suppressed relative to thermal noise. The map 
generated by this technique, illustrated in Fig. [14] shows 
that for the PGB-8 portion of the map, noise levels range 
from 10 mJy to 50 mJy (620 mK to 3.1 K using a syn- 
thesized beam area of 2.15e-5 steradians at 156.4 MHz) 
in the declination range of peak sensitivity of the array. 

6.2. Thermal Noise Contributions to PGB-8 Map 

The map statistics in the direction of a colder patch 
of the galactic synchrotron were evaluated using 1.46- 
MHz facet images with phase centers at (11^*20, 30°00), 
(ll''20,40°00), (12''00,30°00), and (12^00, 40°00). Each 
frequency interval and phase center was imaged twice 
with independent sets of alternating (odd and even time 
index) 14-second integrations from the 3-day PGB-8 ob- 
servation described in |6.1[ For each facet, these two im- 
ages were summed and differenced to produce images of 
the celestial sky and of thermal noise, respectively. For 
the following analysis, image statistics were evaluated 
independently over 10° diameter circular patches around 
each phase center and averaged over facets afterward. 

Thermal noise levels were inferred from measured RMS 
pixel values according to the equation: 



2 fciJ 71-013 



(5) 

A2 / X^^Ai^tN{N -1)' 
where fig is the solid angle of a synthesized beam; il^ is 
the solid angle of the primary beam; Av is the bandwidth 
of data used in the image; r is the integration time; N is 
the number of antennas in the array; and < ... > denotes 
the RMS value. The thermal noise images produced by 
differencing integrations yield the Tims estimates plotted 
in Fig. 11 for each 1.46-MHz band. 

From each measured Trms we also inferred the system 
temperature Tgys present before integration in the corre- 
lator. A comparison of this temperature with the Tgky 
predicted for a colder patch of the galactic synchrotron 
(see [3.2 \ shows that the PGB-8 system temperature is 
dominated by sky-noise. Finally, the standard deviation 
of pixel values in the summed facet images exceeds ther- 
mal noise, implying contributions from celestial point 
sources, their sidelobes, and galactic synchrotron emis- 
sion. The temperature associated with this pixel dis- 
tribution reflects a sky temperature as perceived by the 
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Fig. 11. — We imaged a colder patch of synchrotron sky using 
3 days of PGB-8 data grouped into 1.46-MHz bands, adding inte- 
grations with alternating signs so that thermal noise could be esti- 
mated from RMS pixel values (solid, triangles). From these values 
a system temperature was inferred (solid, dots) using Equation 
Comparison with a predicted synchrotron spectrum for this pati 
of sky (dot-dashed) indicates that sky-noise dominates the system 
temperature of PGB-8. The average sky temperature (dotted, 2cr 
error bars) measured over these facets is substantially lower than 
the predicted synchrotron spectrum because large-scale structure 
is resolved out by the partially filled aperture of PGB-8. 

partially filled aperture of PGB-8. This perceived tem- 
perature is approximately a factor of 6 lower than the 
temperature of the galactic synchrotron because sparse 
samplings of short MiJ-spacings cause much of the emis- 
sion from this source to be resolved out. 

6.3. Angular Power Spectra 

The angular power spectra illustrated in Figure 12] are 
measured using the PGB-8 data described in §6 Our 
analys is follows the techniques developed in [White et al.| 
( 1999 1 for calculating angular power spectra with inter- 



lerometric data in the flat-field approximation. To make 
this approximation, gridded visibility data projected to- 
wards each of the four phase centers defined in |6.2| were 
limited to a 10° diameter FoV by applying a circular box- 
car windowing function W, whose effect in the ut;-plane is 
to convolve by the complementary Fourier kernel Wi^il). 
Within each facet, angular power spectra were calculated 
using the equation: 



C{u) 



2i 



w 



2t:u 



(6) 



where fi-yv represents the solid angle spanned by W; 
and V{u) represents the measured visibility as a func- 
tion of position u = {u, v) in the uu-plane; f denotes 
integration around a ring of |u| — u. The relationship 
between ui;-coordinates and angular wavenumber is de- 
fined by £ ~ 27r l-uj. In practice, V{u) is accompanied by 
sample weights B{u), so that Equation [6] becomes: 
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Fig. 12. — Shown above are angular power spectra measured at 
146.9, 155.7, 164.5, and 173.3 MHz (solid, top to bottom at £ = 100, 
respectively, with 2a error bars), averaged over 1.46-MHz intervals 
for four 10° facets near ll''40, 35°00. These spectra contain contri- 
butions from both point sources and galactic synchrotron emission. 
Dominant errors are from cosmic variance, which has been esti- 
mated from the standard deviation of measurements at the same £ 
in the uv-plane. Dot-dashed lines indicate angular power spectra of 
"noise images", generated by adding integrations with alternating 
signs, and are indicative of the thermal noise level in these mea- 
surements. Also plotted: a fiducial reionization power spectrum 
at z = 9.2 {u =140 MHz, thick solid) and a predicted galactic 
synchrotron spectr um at 147 MHz (dotted) ( [Tegmark et al.|2000| 
[Santos et al.|2005| |. 

where fl^ represents the solid angle spanned by the pri- 
mary beam. 

For each facet, separate C(u) spectra were formed in 
1.46-MHz intervals by integrating around rings in the 
uv-plane that were logarithmically spaced in radius and 
width. Within each ring, data were weighted by the 
square of the primary beam response toward phase cen- 
ter for optimal SNR. We estimated the contribution of 
thermal noise to the error in each C{u) measurement by 
generating power spectra for thermal noise images. This 
was done using the technique of adding integrations with 
alternating signs described in the previous section. The 
error from cosmic variance was much larger than the er- 
ror from thermal noise, owing to the limited aperture 
sampling of the PGB-8 array and the limited area on the 
sky employed for this analysis. Error contributions from 
cosmic variance were estimated by measuring the stan- 
dard deviation of all C(u) measurements in rings of con- 
stant u in the uf -plane for each facet, and dividing this 
by the square root of the number of independent C(u) 
measurements made. We estimated independent mea- 
surements by calculating the size of a coherence patch in 
the uv-plane that arises from W(u) and then degrading 
the pixel resolution in the uv-plane to a corresponding 
size. At the lower resolution, each pixel sample was con- 
sidered an independent measurement, with the caveat 
that for a real- valued sky, only half of the samples in a 
|m| = u ring are independent. 

PGB-8 angular power spectra are dominated by un- 
clustered point sources, as indicated by the approximate 

power-law in £{£ -\- l )Ci /2tt of the angular power spec- 
tra plotted in Figure |12[ However, at lower radio fre- 
quencies for £ < 100, we see a flattening of these spectra 
to a shallow, negative power-law that shows qualitative 
agreement with the C^^sync oc * model adopted by 
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Tegmark et aL] ( |2000[ ) and [Santos et al\ ( |2005[ ) . Progress 
in modeling point sources and in subtracting spectrally 
smooth emission associated with synchrotron sources will 
facilitate the suppression of these foregrounds, with a 
fundamental limit pla ced by the thermal limit of these 
observations (Fig. 12 dot-dashed) , which we measure to 
be 310 mK for £ =^00 at 160 MHz. Overall, our mea- 
surements indicate that over the range 10^ < £ < 10^, 
the fundamental sensitivity level of these observations is 
between two and three orders of magnitude in tempera- 
ture above the fiducial model of reionizatio n at z = 9.2 
{ly =140 MHz) used by [Santos et~al] (|2005 1. 



7. CONCLUSION 

PAPER is following an incremental build-out strategy 
that emphasizes the characterization and optimization of 
instrumental performance at each stage of development. 
We have given special attention to designing antenna el- 
ements whose smooth spatial and spectral responses are 
conducive to calibration that will be required to model 
and remove strong foregrounds to the 21cm EoR signal. 
Facilitating the incremental build-out of PAPER is a se- 
ries of digital correlators that take advantage of the scala- 
bility and reusability of CASPER's packet-switching cor- 
relator architecture. We are developing post-correlation 
calibration, imaging, and analysis pipelines as part of the 
open-source AIPY software package. This package em- 
phasizes the needs of low-frequency interferometry and 
facilitates experimentation with new analysis architec- 
tures and calibration techniques. 

With results from the PWA-4 and PGB-8 deploy- 
ments in Western Australia and Green Bank, respec- 
tively, we demonstrate a first level of calibration that 
is sufficient for modeling strong sources to a level where 
ionospheric refraction and temperature-dependent gains 
must be taken into account. This calibration relies heav- 
ily on a static model of the primary beam. Next steps 
in calibration will focus on the first-order removal of 
temperature-dependence in antenna gains and on mea- 
suring deviations of the primary beam of each antenna 
from the computed model. 



Data from the PWA-4 and PGB-8 deployments have 
been used to generate an all-sky map that attains a ther- 
mal noise level of 10 mJy/Beam (corresponding to 620 
mK, for a 2.15e-5 steradian synthesized beam size at 
154.4 MHz) integrated across a 138.8-MHz to 174.0-MHz 
band. This achievement represents a first iteration in a 
cycle of improvement wherein sky models are used to 
improve array calibration, which in turn allows increas- 
ingly accurate sky models to be generated. In analyzing 
the noise characteristics of this map, we find the sys- 
tem temperature of the PGB-8 array to be consistent 
with a model of galactic synchrotron emission. Angular 
power spectra generated from these data indicate that 
point sources are currently the dominant foreground to 
the EoR signature, with evidence for contributions from 
galactic synchrotron emission at lower radio frequencies 
and angular wavemodes. Errors in measurements are 
dominated by cosmic variance, but by differencing inte- 
grations, we measure the thermal noise level at £ = 100 
and 160 MHz to be 310 mK-a sensitivity threshold that 
lies between two and three orders of magnitude in tem- 
perature above fiducial reionization models. These re- 
sults demonstrate the need for a next level of calibra- 
tion, modeling, and foreground suppression that will be 
pursued in the next PAPER deployments. 
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